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ABSTRACT 



The (n, 7)-cross sections for neutron-rich oxygen isotopes have been calculated in the direct capture model. The 
experimental data for ^^0(n, 7)"'^^0 can be reproduced using this model. Compared to previous work the cross section 
for ^^~^^0(n, 7)^°~^^0 are enhanced considerably by factors between four and 300. 



1. INTRODUCTION 



In the last years the importance of the direct-reaction (DI) mechanism in nucleosynthesis has been realized. 
The DI may even dominate over the compound-nucleus (CN) reaction mechanism in nuclear processes relevant for 
nucleosynthesis. This can be the case for primordial and stellar nucleosynthesis, f. i. in neutron-induced reactions 
in the inhomogenious big-bang scenario^^'^, in many solar nuclear reactions*, in the famous triple-alpha process 
leading to the creation of carbon^, or in miclear reactions of the r-process involving magic-shell nuclei and nuclei far 
from stability^. 

In this work we calculate cross sections of neutron-capture on oxygen isotopes. The reaction rates for these 
reactions are of importance in the inhomogeneous big-bang scenario^'^ as well as in the a-rich freeze out of the 
neutron-rich hot neutrino bubble in supernovae type 11^'^". We investigate the capture reaction "'^^0(n,7)"'^^0 at 
thermal (0.025 eV) and thermonuclear (30keV) projectile energies and compare the calculated cross sections with 
the experimental data. Using the same methods we also calculate the cross sections for the neutron-capture reactions 
i9-2iQj-j-^^^^20-22Q thermal and thermonuclear energy range. 

In section 2 the folding procedure for the determination of the optical and bound-state potentials is introduced. 
In section 3 we discuss the model necessary for the calculation of direct capture (DC) cross sections. The method of 
bound states embedded into continuum (BSEC) allows including isolated CN-resonances in the DC-cross section. 
In section 4 the calculations and results for the neutron capture cross sections for some neutron-rich oxygen isotopes 
are presented and discussed. Finally, in the last section a short summary is given. 



2. FOLDING PROCEDURE 



The folding procedure is used for calculating the nucleon-nucleus potentials in order to describe the elastic 
scattering data and the bound states. This method was already applied successfully in describing many nucleon- 
nucleus systems. In the folding approach the nuclear density pA is folded with an energy and density dependent NN 
interaction VcK^^'^"^ 

V{R) - \Vy{R) ^xJ pAir^vME,PA, \R'r^)dr (1) 

with R being the separation of the centers of mass of the two colliding nuclei. The normalization factor A is adjusted 
to elastic scattering data and to bound- and resonant-state energies. The potential obtained in this way ensures the 



correct behavior of the wave functions in the nuclear exterior. At the low energies considered in nucleosynthesis the 
imaginary parts of the optical potentials are small and can often be neglected. 



In the folding approach the nuclear densities pA for stable nuclei are derived from experimental charge distribu- 
tions^^. For unstable nuclei these densities are determined with the help of the relativistic mean field theory (RMFT). 
The RMFT describes the nucleus as a system of Dirac nucleons interacting via various meson fields. In the last few 
years this theory has turned out to be a very successful tool for the description of many nuclear properties (for 
example binding energies and charge radii for stable isotopes)^'^*. 



In nuclear reactions two extreme types of reaction mechanisms can exist: the compound-nucleus (CN) and the 
direct (DI) process. In the CN mechanism the projectile merges in the target nucleus and excites many degrees of 
freedom of the CN. The excitation proceeds via a multistep process and therefore has a reaction time typically of 
the order 10^^® s to 10"^" s. After this time the CN decays into various exit channels. The relative importance of 
the decay channels is determined by the branching ratios to the final states. In the DI process the projectile excites 
only a few degrees of freedom (e.g. single-particle or collective). The excitation proceeds in one single step and has 
a characteristic time scale of 10"^"'^ s to 10~^^s. This corresponds to the time the projectile needs to pass through 
the target nucleus; this time is much shorter than the reaction time of CN processes. 

In thermonuclear scenarios the projectile energy is well below the Coulomb and/or centrifugal barrier. At these 
energies the competition between different reaction mechanisms is quite complicated. At these energies the CN 
formation may be suppressed, because there may exist no CN levels that can be populated, especially in light, magic 
and far-off-stability nuclei. 

In this work the theoretical cross sections are given through the DC-contribution. The theoretical cross section 
cr*'' is obtained from the DC cross section a^'~^ by^'^^ 



The sum extends over the ground state and excited states in the final nuclei, where the spectroscopic factors Si 
are known. The isospin Clebsch-Gordan coefficients are given by d. The DC cross sections af"~^ are essentially 
determined by the overlap of the scattering wave function in the entrance channel, the bound-state wave function 
in the exit channel and the multipole transition-operator. The radial dependence of the DC-integral is in our case 
determined uniquely by the folding potentials. 

Isolated CN-resonances with known resonance energy and width can be incorporated into the nonresonant DI 
cross section by BSEC^®""'^^ which allows a simultaneous calculation of the resonant and non-resonant contributions 
of the cross sections. In this approach the scattering wave function is calculated using an adequate energy-dependent 
potential generating single-particle resonances and reproducing effects caused by their coupling with complicated 
bound states into the continuum. 



The density distributions of ^^0 to ^^O were derived with the help of RMFT.^'* In the entrance channel the 
normalisation factor A given in Eq. |l| of ^^0+n was adjusted to reproduce the total cross section. For the other 



3. REACTION MECHANISMS AND MODELS 




(2) 



4. CALCULATIONS AND RESULTS 



Figure 1: Level scheme for the reaction ^*0(n,7)^^0 and comparison of the DC cross section for ^*0(n,7)^^0 with 
the experimental data^'*'^^'^^. 

reactions the A were adjusted to reproduce the same volume integral J = 578.33 MeV fm'^ as for the folding potential 
for ^^0+n elastic scattering. In the exit channels the A were fixed by the neutron separation energy from the final 
nuclei. The folding potentials of Eq. |l| were determined with the help of the computer code DFOLD.^" 

The spectroscopic factors for one-nucleon transfer of the oxygen isotopes were determined from shell-model cal- 
culations in the sd-shell with an effective nucleon-nucleon interaction derived by Wildenthal^^. For these calculations 
the programs GENESIS, RITSSCHIL, DIA22 and SPECTROS^^ were used. The resonance energies and widths of 
^^O used for the BSEC-calculations of ^^0(n,7)^^0 are listed in Ref. 24. For the direct-capture calculations Eq. || 
the code TEDCA^^ was used. For the BSEC-calculations the code TEDCA was generalized to include this approach. 

Let us first discuss the reaction ^®0(n,7)^^0. For this reaction experimental reaction data at the thermal energy 
(0.025 meV) and around the 3/2~-resonance at 625 keV are available. The value of the cross section at 23.3keV is 
given by ~ 8 ± 1 /ibarn^^'^^. We consider this reaction as a test for our DC-model. 

The level scheme for this reaction is shown on the left side in Fig. 1. There are two types of El-transitions 
relevant for the energy range between thermal and thermonuclear energies. The first one is from an s-wave in the 
entrance channel exciting the negative-parity state 3/2^ in the final nucleus which is just bound by about 12.1 keV. 
Even so this transition has a very low Q-value it still reproduces the experimental cross section (Jthcrm = 0.16mb at 
0.025 eV with a reasonable spectroscopic factor of 0.331. In this case the transition matrix has contributions up to 
500 fm and has a maximum at about 100 fm. This transition gives the well-known 1/v-behavior (see Fig. 1). The 
second type of transition comes from an initial p-wave and excites the positive-parity states 5/2"'', 3/2+, 1/2"'" in 
the final nucleus. These transitions have a v-behavior (see Fig. 1). Furthermore as can be seen from Fig. 1, the 
resonance 3/2"" at 625 keV can be described by the BSEC which was discussed before. 

The contributions of the different transitions to the cross section at 30keV and their sum for the ^^~^-'^0(n,7)-'^^~^^0 
reactions are shown in table 1. Also the spectroscopic factors obtained from the shell-model calculations are listed. 
One can see, that the cross sections calculated in this work are enhancend by factors between four and 300. The 
reasons are that in previous work® spectroscopic factors have been assumed too low and important transitions to 



Table 1: Comparison of the contributions of different transitions to bound-states and their spectroscopic factors 
between previous work^ and this work for the reactions ^^~^^0(n,7)^^~^^0 



reaction 


final state 


spectroscopic factor 


transition 


previous work 


this work 












Sn [/Ltb/VMeV] {E 


= 30keV) 


i*0(n, 7)1^0 


5/2+ 


0.687 


P - 


■^d 




8.85 




3/2+ 


0.013 


P - 


^d 




0.13 




1/2+ 


0.830 


P - 


-> s 




44.43 








total 


52.40^ 


53.33 


^y0(n,7)^"0 


0+ 


3.427 


P - 


■^d 


0.54 


2.72 




2+ 


0.731 


P- 


d 


1.94 


2.63 






0.142 


P - 


s 


4.15 


5.16 




4+ 


1.021 


P- 


^ d 




4.63 




2+ 


0.574 


P - 


s 




15.58 








total 


6.63 


30.72 


20O(n, 7)210 


5/2+ 


0.345 


P - 


^d 


11.02 


5.40 




1/2+ 


0.811 


P - 


-> s 




40.23 








total 


11.02 


45.63 


^iO(n,7)^^0 


0+ 


5.222 


P - 


■^d 


0.16 


4.81 




2+ 


0.821 


P - 


s 




22.42 




3+ 


0.771 


P - 


s 




21.44 








total 


0.16 


48.67 



^ Determined with the help of the experimental value of the cross section at 23.3 keV. 
some excited final states have been neglected. 



5. SUMMARY 



Direct-capture calculations using the folding procedure can reproduce excellently the experimental data for 
neutron-rich nuclei in the thermal as well as in thermonuclear energy region. Isolated resonances can also be 
incorporated sucessfully in the direct capture model. The cross sections calculated in our model for neutron -rich 
oxygen isotopes are for some reactions considerably higher than the ones given previously. Such changes should lead 
to modifications of existing reaction network calculations. 
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